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Abstract 

Data collected by the KLOE experiment in years 2001/2002 for a total integrated 

luminosity of 410 pb~^, corresponding to 17 million of produced r/'s from the $ ^ 777 
radiative decay, have been analyzed to look for the C- violating decay rj — > 777. The 
signal is searched in events with 4 photons in the final state by using the spectrum 
of the most energetic photon. The background is evaluated from the data themself, 
doing a fifth degree polynomial fit outside the signal region, and extrapolating 
inside it. No signal has been observed and taking into account a selection efficiency 
e = 20.3% we get from a likelihood fit Br(?? 777) < 1.6 x 10"^ at 95% CL. The 
procedure has a systematic ambiguity at the 35 % level. 



1 Introduction 

At the Frascati factory DA<I>NE [1] a large number of rj mesons is produced 
via the radiative decay of the 4> meson. 
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The mesons are produced in the colhsion of e+e" beams into the DA$NE 
colhder and 1.298 % of them decay into the rij final state. Up to now 17 x 10^ 
rj mesons have been produced and are present in 410 pb~^ collected by the 
KLOE experiment in years 2001-2002. 

The 1] meson is an even eigenstate of the charge coniugation operator C, 
while a photon is an odd C eigenstate. So any decay of the rj meson into 
a final state with an odd number of photons violates the C symmetry. In 
the Standard Model, the C symmetry is exactly conserved in both strong 
and electromagnetic decays, but it is violated in weak decays due to the V-A 
structure of the weak couplings. In the framework of the Standard Model the 
decay rate of the 7r° — > 87 has been evaluated [2], and generalizing this result 
to the rj case [3], one obtains Br(77 — > 87) < 10~^^. 

For this reason the discovery of a larger decay rate would be a clear signal of 
SM violation. At the moment all the predictions of alternative models are far 
below the experimentally achievable limits [3] . 

Prom the experimental point of view the only published result is that of 
GAMS2000 experiment (1984) which has obtained the upper limit 5 x 10"'^ 
at 90 % CL. [4]. There is also a preliminary result from Crystall Ball collab- 
oration (AGS/CB) at Brookhaven which sets the Umit at 1.8 x 10~^ at 90 % 
C.L [3]. Both the above experiments are of fixed target type, in which the rj 
meson is produced through the reaction Tr~p — > rjn with P^- — 30 GeV/c for 
GAMS and 760 MeV/c for AGS/GB. 



2 Kloe detector 

The detector consists of a large cylindrical drift chamber, DC [5], surrounded 
by a lead-scintillating fiber sampling calorimeter, EMG [6], both immersed 
in a solenoidal magnetic field of 0.52 T with the axis parallel to the beams. 
The DC tracking volume extends from 28.5 to 190.5 cm in radius and is 340 
cm in length. For charged particles the transverse momentum resolution is 
6pt/pt — 0.4% and vertices are reconstructed with a spatial resolution of ~ 

3 mm. The calorimeter is divided into a barrel and two endcaps and covers 
98% of the sohd angle. Photon energies and arrival times are measured with 
resolutions as/E = 0.057/^Je (GeV) and at = 54 ps/^E (GeV) 50 ps 
respectively. The photon entry points are determined with an accuracy ai ~ 
1 cm/ ^E (GeV) along the fibers, and ~ 1 cm in the transverse directions. A 
photon is defined as a calorimeter cluster not associated to a charged particle, 
by requiring that the distance along the fibers between the cluster centroid 
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and the impact point of the nearest extrapolated track be greater than 3ai. 
Two small calorimeters, QCAL [7], made with lead and scintillating tiles are 
wrapped around the low-beta quadrupoles to complete the hermeticity. 

The trigger [8] uses information from both the calorimeter and the drift cham- 
ber. The EMC trigger requires two local energy deposits above threshold 
{E > 50 MeV in the barrel, > 150 MeV in the endcaps). Recognition 
and rejection of cosmic-ray events is also performed at the trigger level, check- 
ing for the presence of two energy deposits above 30 MeV in the outermost 
calorimeter planes. The DC trigger is based on the multiplicity and topology 
of the hits in the drift cells. The trigger has a large time spread with respect 
to the beam crossing time. It is however synchronized with the machine radio 
frequency divided by four, Tgync = 10.85 ns, with an accuracy of 50 ps. Dur- 
ing the period of data taking the bunch crossing period at DA$NE was T — 
5.43 ns. The To of the bunch crossing producing an event is determined offline 
during the event reconstruction. 



3 Background 

The signal in this analysis is — > 7?y(— >■ 87) — >■ 47, the cross section corre- 
sponding to the PDG2002 upper limit is very low: < 0.0215 nb (at KLOE 
the 4> production cross section is about 3 /^b). So, background studies should 
cover all possible neutral processes at KLOE. The most important process 
which gives 4 photons in the final state is e+e~ — > cu^ with cu — > 7r°7 but also 
important are processes with less or more than 4 photons in final state which 
mimic 4 photons events because of cluster splitting or merging and accidental 
clusters due to machine background. The agreement with MC in the back- 
ground description is very hard to obtain due to the above effects which are 
very difficult to reproduce; so we have estimated the background directly from 
data and used the MC only to evaluate the detection efficiency of the signal. 

To this aim a generator of 77 — >777 with pure phase space for the decay 
dynamics has been used to produce 120,000 — > 777(77 — > 37) events. 



4 Event selection 



First a preselection is applied to the reconstructed events of the neutral ra- 
diative stream. A 'recover splitting' procedure (applied both to data and to 
MC events) is implemented to reduce the number of split clusters; moreover 
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it is rcqTiircd that every cluster does not have an association with a DC track. 
Reconstructed data also include beam position [8] and energy momentum 
determination, obtained run by run from the analysis of Bhabha events. MC 
includes cluster time offset correction, and simulation of software filters used 
to select events for the neutral radiative stream. 

The following cuts are applied to the events: 

1) Ntime=4 

Requires 4 photons with the correct time of flight t : \t—r/c\ < 5(Jt ; r being the 
distance of the cluster from LP. and at being the calorimeter time resolution. 

2) Ngood = 4 

A good photon is defined in the following way: cluster energy > 30MeV, 
\cos6\ < 0.93, where 6 is the angle respect to the beam direction. Cluster 
energy cut helps to reject accidental clusters and split clusters, while the polar 
angle cut excludes the blind region around the beam-pipe. 

3) Total energy of prompt clusters > 800 MeV 

4) Total momentum < 200MeV 

These last two cuts reject events from channels with more than 4 photons in 
the final state which give background if some photons are lost. 

The events that pass the above preselection criteria are stored in reduced files 
and to them are applied the following additional cuts: 

5) e^j>i5° 

9^^ is the minimum angle between two photons, to further reject split clusters, 
which mainly come from 37 final states. 

6) E^i^ > 50 MeV, \cose\ < 0.91. 

These cuts further reject accidental clusters and QED background. 
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Entries 61^?7 

Meon 10.54 

RMS 6.906 



20 30 40 50 60 70 80 90 100 



Fig. 1. of kinematic fit distribution for: top DATA, bottom MC signal. Events 
with > 25 arc rejected. 

5 Kinematic fit 



At this point a kinematic fit procedure is applied to the four photons to 
improve the energy resolution. The input variables of the fit are the following: 

• X,Y,Z coordinates of the cluster; 

• E energy of the cluster; 

• t time of the cluster; 

• X,Y,Z of the interaction vertex; 

• E,Px,Py,Pz of the meson. 



The fit is made according to the Lagrange's multipliers method, minimizing 
the following x^: 



X 



with 27 free parameters, while the Fj[iik) represent the energy, momentum 
and time constraints. 



The of the kinematic fit is used to reject the background using the cut 
X^j„ > 25 (see Fig. 1). At this point the main source of background is given 
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by the channel e+e^ uj{^ n^j)j where a photon comes from ISR radiation. 
This is evident by making a kinematic fit of the event in the 71*^77 hypothesis 
and looking to the plot of the m^o^ variable (choosing the most energetic one 
between the two photons not linked to the 7r°) (see Fig. 2). 




400 600 




1 uuu 
MeV/c' 



Fig. 2. 



distribution, it shows a clear peak at the to mass. 



Events with a tt" in the final state represent a very large fraction of background, 
they come also from /07 7r''7r°7, from (j) 0,07 rjn^'-f and —> 7r°7 
final states. We reject these events containing a 7r° by cutting on the invariant 
mass m^^ built from any couple of photons. This variable is plotted in Fig. 3 
both for 2001 data and MC signal. The cut chosen is 90 MeV/c^ < m^^ < 
180 MeV/c2. 



6 Upper limit evaluation 



To search for the r] 87 events, we look for a peak in the distribution of the 
energy Emax (evaluated in the reference frame) of the most energetic photon 
. For the majority of the events this photon corresponds to the radiative one 
which has an expected energy of about 363 MeV. In Fig.4 is shown the Emax 
distribution for MC, 2001 and 2002 data. 

As can be seen the two distributions for 2001 and 2002 overlap very well. 
A Kolmogorov test has been performed to check their compatibility, it gives 
a compatibility probability of 26%, so we use the whole sample. From the 
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Fig. 3. rrir^r^ distributions. Top: 2001 DATA, bottom: MC signal. Events in the 90 - 
180 MeV mass window are rejected. 




Fig. 4. Spectrum of the most energetic photon in the (j) reference frame: left: 2002 
data (red histogram) and 2001 data (blue points); right: MC signal r/ ^ 87 

distribution is very clear that there isn't a narrow peak in Emax, so we don't 
see any evidence of 77 — 87 events. 



Then to evaluate the upper limit on the — >^ 87 signal using this distribution 
we choose as signal region (the region where there is the main part of the sig- 
nal) the interval [850,880] MeV. Then we fit the observed distribution of E^ax 
in the domain [280,850] MeV U [880,480] MeV with a fifth degree polynomial. 
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The fit is good (x^/n.d.o.f = 77.9/92); (see Fig. 5) and the fitted polynomial 
is used to estimate the number of the background events in the signal region 
[350,380]. 




300 320 340 360 380 400 420 440 460 480 
E,.„, MeV 




300 320 340 360 380 400 420 440 460 480 
E„„ MeV 



Fig. 5. Top: Emax distribution for 2001+2002 data. The interval used for background 
estimation is shown; bottom: the fitted polynomial is overlapped to the whole range. 

In this range we assume to have both background and signal and build a 
likelihood function in this way: 



Likelihood{N, 



sigj 



n 

i=l 



1 



27rN: 



exp 



■ exp — 



2Nr' 



where Ni is the number of observed events in the bin "i" (we use in the signal 
region 15 bins of 1.75 MeV width), and N^^^ is the number of the expected 
events given by: 



N 



exp 



ai + a2 ■ Xi + as ■ x^^ + ■ + ■ xf + ttQ ■ xf + Ngig ■ f, 



where Xi is the central value of the bin "i" and fi is the fraction of signal events 
that fah in the bin "i" /j = !)• The /« are evaluated using the histogram 
of Fig. 4(right). The likelihood function and its integral function are shown in 
Fig. 6. We obtain the following values for the upper limit on Ngig: 

Nsig < 59.0 at 95 % C.L. 

Ns^a < 49.2 at 90 % C.L. 



9 




-26 



40 50 60 70 80 90 100 



1 

0.8 
0.6 
0.4 
0.2 




- 0.9b 





















10 20 30 40 50 60 70 80 90 100 



Fig. 6. Top: Likelihood as a function of Ngig] bottom: Integral function of the like- 
lihood; the 95 % and 90 % limits are shown. 

To convert the above results into the corresponding upper limits on the Br (77 —>■ 87) 
we have evaluated the number of r] in the data sample through the study of the 
decay channel rj — > Svr*^ by using an already developed analysis inside the ex- 
periment [9] . The efficiency of the selection evaluated from Monte Carlo is e = 
20.3 %. The result is: iV(r/ ^ 371°) = 2431917 and e^^g^o = 0.4106 ± 0.0016. 
Then we evaluate the ratios of the two branching ratios: 

Br{r] 37) ^ Nup e^^g^o 



77^37 



Br{r] — > 37r") 



Br{rj —>■ 37) 
Br{ri 37rO) 



< 4.1 ■ 10" 



90 % C.L. 



Using the PDG2002 value for Br{r] 37r°) = (32.51 ± 0.29)% we have: 
Br{ri 37) < 1.6 • 10"^ 95 % C.L. 
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Br{rj 37) < 1.3 • 10"^ 90 % C.L. 



7 Systematic errors 

In this section we evaluate various systematic effects, especially those which 
may arise from differences between the data and MC distributions used in the 
analysis. 



7.1 Systematics due to the x^cut 

We evaluate the systematics due to possible differences in the distribution 
between data and MC by comparing a data control sample with MC. The 
chosen control sample is the channel e^e^ — a;(— > tt^'j) — >■ 47 that is the only 
channel with four photons that we have in our selection. Since the kinematic 
fit requires only the energy-momentum conservation, we can compare its 
distribution directly to that of the channel rj ^ 3'j from MC. We select the 
uj'-f candidates by requiring the mass range 128 < m^^ < 145 MeV/c^ for 
the 7r° and 760 < m^Oy < 815 MeV/c^ for the u. In Fig. 7 we have reported 
the distribution of the kinematic fit, under energy- momentum conservation 
hypothesis, for the selected cu^ sample and for the MC 77 — > 37 sample. The 
upper plot is the x^ distribution, the lower plot is the fraction of events that 
survives to a given x^ cut, normalized to the range shown in Fig. For xlut — 25 
the data-MC discrepancy is about 3%. This is the systematic error that we 
assume for this cut. 



7.2 Correction for photon detection efficiency 

The MC doesn't simulate with full accuracy the photon detection efficiency, 
from the study done by [10] we know, as a function of the energy of the photon 
and the cos{9), the ratio w = edata/^MC- For this reason, we have evaluated 
the quantity: 

4 

Weight = Wi "i" runs on the four photons. 

i=l 



11 
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Mean 10.47 

RMS 6.570 




Fig. 7. Top: distribution for MC rj 87, continuous line, and uj selected sample 
(see text), dashed line. Bottom: integral of the above distributions. 

for each event. The sum of the weights gives the effect of this discrepancy on 
the efficiency. The results are: 



unweighted events 


24376/120000 (20.3 %) 


weighted events 


24167/120000 (20.1 %) 


Ae/e 


1.0 % 



7.3 Systematics due to signal or background shape and window choice 



To test if the kinematic fit introduces a bias in the photon energy that is 
different for data and MC, we have analyzed a sample of 37 events (0 

77)7) and obtained from the kinematic fit the updated energy for every 
photon. This energy is plotted in Fig. 8 both for data and MC. After having 
taken in account the wrong rj mass that we have in MC generator (m^ = 548.8 
MeV instead of 547.3 MeV) the two distribution overlap very well, so we don't 
quote the systematics on this variable. 

To take into account other possible systematic effects in signal and background 
shape we varied the width of the signal window and the order of the polynomial 
used to fit the background shape. Then we repeated the fit of the background 
(background+signal) excluding (including) the new energy window and by 
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Radiative piioton energy residual MeV 



Fig. 8. E^ — E^rad distribution for photons in the (f) — > 77)7 sample. Each event 
enters three times in the histogram, one for each photon. Both MC (continuous 
hne) and data (points) are shown. E^rad is the expected energy of the radiative 
photon = 362.7 MeV. The MC shape has been fitted with a double gaussian fit: 
PI ■ exp[-^:^] + P4 • exp[- (^#]. 

recalculating also the signal efficiency: the maximum variation found on the 
various upper limit estimates was about 35 %. 



8 Conclusions 

We have obtained the upper limit Br{ri — * 37) < 1.6 ■ 10~^ at 95 % C.L. 
which is significantly better than the PDG existing one and also improves the 
unpublished Crystal Ball value. We estimate our systematic uncertainties to 
be at the 35% level. More data to be collected by Kloe in the next future will 
allow to further improve this result. 
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